Cognitive Neurodynamics

Abnormal states such as erroneous behaviors are generally difficult to represent from neural data. However, such states are also known to have specific spatiotemporal features, indicating a feasibility of developing a method to focus on them. If a method can highlight these spatiotemporal features, it may effectively represent such abnormal states, helping evaluate abnormal brain functions. In the present study, we proposed the hierarchy of supported modules (HSM) to highlight spatiotemporal features that can represent abnormal states. HSM spatiotemporally transforms multidimensional neural time-series based on their spatiotemporal context. We evaluated HSM through decoding and similarity analyses using multiple publicly available datasets. In the HSM results, decoding accuracies were higher for erroneous behaviors than for normal behaviors, and similarities were lower between erroneous behaviors and normal behaviors than between normal behaviors, demonstrating the ability of HSM to capture the spatiotemporal features of erroneous behaviors. Surprisingly, many parts of these results were also present even before HSM learning, showing the virtue of HSM as a simple-to-use method. The proposed HSM method may help elucidate the mechanisms underlying erroneous behaviors.

Wave recurrent neural networks (wRNNs) are biologically inspired recurrent architectures that use traveling-wave dynamics to support sequence learning and memory. However, their input-to-hidden pathway remains relatively simple compared with biological neurons, where dendrites perform nonlinear input integration. In this study, we introduce the Dendritic Wave Recurrent Neural Network (DWRNN), which augments the input pathway of the wRNN with nonlinear basal dendritic branches while preserving the original recurrent wave dynamics. We evaluate DW-RNN on a simple copy task, sequential MNIST (sMNIST), permuted sequential MNIST (psMNIST), and noisy sequential CIFAR-10 (nsCIFAR-10). On the copy task, DW-RNN shows learning behavior comparable to the standard wRNN, suggesting that dendritic input integration does not disrupt the recurrent wave-based memory mechanism. On the three sequential image-classification benchmarks, DW-RNN outperforms the standard wRNN, improving accuracy from 97.27 {+/-} 0.15% to 97.82 {+/-} 0.12% on sMNIST, from 96.74 {+/-} 0.17% to 96.92 {+/-} 0.10% on psMNIST, and from 54.30 {+/-} 0.79% to 55.65 {+/-} 0.55% on nsCIFAR-10. In addition to improving mean accuracy, DW-RNN exhibits lower across-seed variability on all three classification benchmarks, suggesting that dendritic input integration may improve the stability of wRNN training. Hidden-activity visualizations further show that DW-RNN preserves the characteristic traveling-wave patterns of the original wRNN. These results suggest that dendritic computation and traveling-wave recurrent dynamics provide complementary mechanisms for biologically inspired sequence learning.

Despite extensive study of cellular mechanisms underlying long-term potentiation, no single specific protein or gene has been identified which encodes an individual unit of information, or memory bit. Indeed, the brain engram remains a knowledge gap. The theory of exclusion led us to cancel one-by-one several unrealistic biological options, suggesting that the explanation resides somewhere else. Superposition of up to concentric 300 myelin layers, spiraled, and highly compacted wrapping a single axon and each wrap could host hundreds to thousands of niches, as memory cells, collectively consisting of a massive array of cells. The disjointed 3D spatial superposition allows storage of charges, nodes not facing from a layer to next. The thickness of a single myelin layer ranges from 7.0 to 20 nm. The dimension scale is approximately the exact dimensions of the charge trap, the tunnel and dielectric also equipping current AI microchips. Stored charges are positive ions, with similar effect whether charges are negative or positive charges creating an electromagnetic field. To write data, following an action potential, this voltage applies to the control gates of the myelin layers producing an ionic charge injection. This causes charges to gain energy and tunnel through the myelin layer across Ranvier nodes, via quantum tunneling, and deep into the concentric myelin multilayers. This is creating an insulated trapping of K+ ions isolated from the system. In a long white matter tract bundle, the near-perfect isolation of millions of axons within compressed myelin wrap-ion channel K+/Na+ systems provides quantum coherence and precision of asynchronous firing property. The injected ionic charges (K+) become physically stuck in traps within the myelin layers. The K+ ions may not move freely, completely trapped after AP ceases. Mirroring a single-bit, single-level-cell, a trapped ionic charge (ions K+) may represent a 1, while an empty cell (absence of K+) represents a 0. The trial-and-error process, with a Bayesian inference which may have also been the core evolution of the learning human brain. Based on selected mathematical equations, we analyzed the general scheme on how deep learning may be embedded in the brain

Learning new information in the presence of distracters and changing conditions requires the ability to adapt. In the brain, this adaptive capability has been linked to dynamic interactions between attention and working memory, which enable the selective filtering of irrelevant input while preserving behaviorally relevant information. Specific neural oscillations have been implicated in this process. Here, we introduce a phenomenological data-driven framework for oscillatory network modeling that learns condition-dependent coupling laws directly from neural recordings and enables inference of condition-dependent directed pathways. We apply our approach to magnetoen-cephalography (MEG) data collected while participants performed a working-memory task with and without distracters. Recall dynamics in the non-distracter condition are first modeled using a linear oscillatory network in which each region of interest is represented by two alpha-band harmonic oscillators. We use universal differential equations (UDE), an extension of neural differential equations, to capture distracter-induced changes in coupling laws. Symbolic regression is then used to interpret the modifications identified by UDE as nonlinear functions, and an additional method is proposed to identify the directed pathway from the newly emerging nonlinear terms in the dynamics of brain regions of interest. Despite inter-subject variability, working memory recall data from all four participants examined under distraction showed the emergence of a pathway from the dorsolateral prefrontal cortex (dlPFC) to the primary visual cortex (V1). This finding is consistent with the established role of the dlPFC in cognitive control and suggests that distracter processing recruits a directed interaction from prefrontal to visual regions. More broadly, our results illustrate that combining linear models whose parameters are learned from the data with universal differential equations augmented by interpretability methods enables the identification of condition-dependent coupling laws, their representation as interpretable mathematical functions, and the discovery of candidate directed pathways underlying adaptive changes in oscillatory networks without requiring strong prior assumptions about the underlying mechanisms.

Biomolecular condensates are widespread cellular self-assembled structures with essential functions. There are suggestions of condensates formed by different proteins being near criticality. However, systematic investigation of the criticality of condensates is absent, and critical exponents defining their universality class have not been found. Here, using long-time simulations, we show that condensates exhibit typical critical phenomena, including scale-free spatiotemporal correlations, critical slowing down, divergence of correlation length and dynamic scaling. From these scaling behaviors, a set of critical exponents is determined. Based on dynamic critical exponent, diverse condensates can be divided into two distinct universality classes, arising from differences in their molecular components and interaction types.

Neuroimaging evidence suggests that a distributed network of brain areas, known as the multiple demand (MD) network, is consistently recruited across a wide range of cognitive tasks. The MD network includes specific areas in prefrontal and parietal cortices. However, the fine-grained organization of this network is poorly understood. Here we aim to comprehensively characterize the functional subdivisions within the MD network using a naturalistic fMRI paradigm. 20 subjects were instructed to perform 14 different tasks on a set of movie stimuli. These tasks were designed to target a wide range of cognitive domains including visual, spatial, categorical, emotional, auditory, linguistic, social, and semantic processes. fMRI data were also collected while subjects passively watched the movies. The MD network was first delineated by localizing cortical areas that were significantly more active in all 14 tasks compared to the passive-viewing condition. The principal component analysis was then applied on task-specific responses of cortical points within the MD network. The first component was correlated with activities for all tasks, and its spatial map revealed the core, highly multimodal subregions of the MD network. The other components showed preferences for a subset of tasks. In particular, the second component revealed a sharp distinction between MD regions that were preferentially active in visual/spatial versus linguistic/semantic tasks. A graph analysis on the entire cortex also showed a large-scale distinction between visual/spatial and linguistic/semantic areas, with MD regions linking the two communities of cortical areas. Overall, our results provide new insights into how the MD network and its fine-grained architecture contribute to the human intelligent behavior.

Working memory (WM) is a core cognitive function that supports goal-directed behavior by temporarily maintaining and manipulating information. One of the most widely used paradigms for investigating WM function is the N-back task, and numerous neuroimaging studies have examined load-dependent neural responses using a variety of analytical approaches. However, most previous studies have focused on low-to-moderate load ranges (primarily 0-3-back), and it remains unclear how whole-brain activity patterns reconfigure across a broader range of WM demands, including conditions approaching capacity limits. In the present study, we investigated behavioral performance and whole-brain activity patterns across an extended N-back task ranging from 0-back to 7-back. Behavioral analyses revealed that discrimination sensitivity (d') decreased nonlinearly with increasing WM load, whereas reaction time (RT) exhibited an inverted-U pattern, peaking at intermediate load conditions. To characterize load-dependent whole-brain activity patterns, we computed relative activation maps by subtracting the participant-wise mean activation map across all conditions from each condition-specific activation map. Spatial similarity analyses with the Yeo 7-network templates revealed that low-load conditions showed relatively high similarity to default mode network (DMN)-related patterns. Similarity to the dorsal attention network (DAN) and frontoparietal network (FPN) was maximal at intermediate load levels, indicating load-dependent changes in network similarity profiles. High-load conditions were characterized by partial re-emergence of DMN-related patterns, accompanied by reduced DAN/FPN similarity. In addition, semantic similarity analysis using Neurosynth-derived semantic maps revealed relatively high similarity to default mode-related and self-referential representations under low-load conditions. Intermediate-load conditions showed strong correspondence with working memory- and executive control-related representations, whereas high-load conditions exhibited increased similarity to salience-, aversive/interoceptive-, and inhibitory-control-related representations. Together, these findings suggest that increasing WM load is associated not merely with stronger activation, but with changes in whole-brain activity patterns accompanied by nonlinear changes in network similarity profiles across levels of cognitive demand. Furthermore, the relative activation map-based whole-brain pattern analysis used in this study may provide a useful approach for evaluating changes in whole-brain state representations associated with cognitive load.

Fast periodic visual stimulation (FPVS) has recently demonstrated an ability to index multiple cognitive domains such as facial expression processing, working memory, and more recently, semantic categorization in just a few minutes of recording time. The present work investigates the effects of low semantic distance and how this modulates semantic categorization responses. Twenty-seven healthy young adults completed an FPVS oddball paradigm in order to determine whether comparing fruits and vegetables of the same color would elicit semantic categorization responses. Strong oddball responses were observed up to the 10th frequency harmonic, and responses showed statistically significant right occipito-temporal lateralization that was not present in a low-level color change condition completed by participants and is opposite the pattern observed in recent word-based semantic categorization FPVS studies. The findings suggest that image-based FPVS paradigms should be investigated further as candidate tools to study conditions that affect semantic categorization such as Alzheimers disease.

Introduction Understanding the mechanisms of human neuromaturation constitutes one of the fundamental questions of neuroscience. While it is well described that large-scale brain maturation is initiated within sensorimotor brain regions and progresses to associative cortex, the underlying developmental neurobiology remains to be fully characterized. Animal models have indicated that cortical inhibitory upregulation might be a driver of neurodevelopment. To investigate the hypothesis that cortical inhibitory upregulation plays a similar role in human neuromaturation, we developed a The Virtual Brain (TVB) based computational model (TVB-Child) to explore potential mechanisms of human neurodevelopment. Material and method We created neurodevelopmental dynamic brain network models capturing neurobiological maturation by using the large-scale brain simulator TVB and fitting brain network models to developmental functional MRI (fMRI) from the Human Connectome Project-Development (HCP-D) data set with 640 subjects with an age range of 6-21 years. Age-dependent trajectories in the fMRI data set were first analyzed by combined group-ICA/Dual Regression extracting subject-specific resting-state networks (RSN). Maturational topographical and topological redistribution of these networks were analyzed by linear and non-linear regression of RSN size and degree and strength centrality. Brain network models were fitted to the fMRI functional connectivity obtained from the HCP-D data set. Hypothesizing that cortical inhibition is a driver of neuromaturation, we analyzed spatiotemporal inhibition parameter gradients in the dynamic brain network model for the hypothesized significant correlations with fMRI RSN maturational trajectories. Results While during development frontoparietal (FP) and default mode network (DMN) grew and exhibited an increase in both degree and strength centrality, becoming dominant network hubs, the attention network underwent network pruning with a decrease in size and node degree. The primary sensory network changed little. For the fitted brain network models, we obtained a high degree of reproduction with correlation coefficients between empirical and simulated functional connectivities ranging between 0.80 and 0.95. Values of the feed forward inhibition model parameter wijFFI representing the strength of regional feedforward inhibitory input exhibited the most significant increase with age within the FP and DMN networks. A less pronounced, but significant, age-dependent increase of the inhibitory parameter values were seen in attention networks and no change within primary sensory networks. Conclusion Our study shows that high order (FP, DMN), attention and primary sensory networks exhibit distinct topographical and topological maturation trajectories. Moreover, brain network modeling revealed RSN-specific age-dependent inhibition trajectories, indicating that the model is able to reproduce and thus support candidate mechanisms of neurodevelopment.

Transformer-based deep learning models have shown great potential for decoding visual EEG signals. However, their internal attention mechanisms are often evaluated primarily on optimization objectives, leaving their alignment with biological brain connectivity an open question. This study empirically evaluates how variations in EEG preprocessing strategies affect these attention representations using the Adaptive Thinking Mapper (ATM) model as a framework. We compared a baseline pipeline (MVNN only) against a comprehensive cleaning pipeline integrating ICA and notch filtering. The models were evaluated through cross-generalization, noise robustness, and spectral-temporal ablation analyses. Furthermore, we investigated the structural correspondence between the model's data-driven attention weights and neurophysiological reference networks (GPDC, PDC, and DTF) using Node Strength Correlation and Representational Similarity Analysis (RSA). The results show that the comprehensive preprocessing successfully suppresses non-neural artifacts, such as frontal noise and electrical interference, while maintaining comparable decoding accuracy and baseline robustness. Alignment analyses revealed that the broad spatial organization of the learned attention patterns remains highly stable across pipelines, capturing key directed connectivity dynamics with subtle, metric-dependent variations in global representational geometry. This work provides an empirical exploration into bridging data-driven attention weights with neurophysiological consistency, offering insights toward more transparent brain-computer interfaces.

This study investigated how different forms of task influence the oral performance of Chinese second-language learners and native speakers. By analyzing data from 40 Chinese second-language learners and 40 native speakers through picture description tasks, formal and informal questions, and questions with different emotions (happy and unhappy), it was found that different task characteristics significantly affected language performance. Short picture tasks led to higher communication efficiency and noun rates but more errors, while long story tasks showed higher verb rates, function word rates, etc. Formal questions had more characters and nouns but lower communication efficiency compared to informal ones. Also, happy emotion questions resulted in fewer characters, sentences, and errors than unhappy emotion questions. These findings contribute to the theoretical understanding of task-based language performance in Chinese as a second language and offer practical implications for teaching, textbook compilation, and student evaluation.

Derived from previous observations on equal and cross-frequency coupling, we evaluated the proposal that equal and cross-frequency phase synchronization may characterize the integration-segregation perspective of cerebral sensory-motor processing. Using brain recordings obtained in normal conditions and in conditions of diminished sensory input (eyes closed wakefulness, sleep and coma, when there is presumably less functional segregation of sensory-motor processing in neural networks), we assessed potential differences in partitioning of the synchrony state space linked to cross-frequency synchronization. More partitions were found in conditions of decreased sensory input. In addition, there was a less complex synchrony state space in cross-frequency as compared with equal-frequency coupling, in terms of fewer connectivity configurations. These results support the idea that equal-frequency coupling favours integration from multiple brain regions occurring in a complex synchrony state space rich in possible connectivity configurations, whereas cross-frequency coupling contributes to segregation, or localized sensory-motor transformations taking place in specific brain areas. This evidence may contribute to new considerations about the much-discussed role of multi-frequency relations in neuronal activity, and how the structural and functional modular organization of the nervous system is able to generate the coordinated activity needed for conscious and appropriate cognitive behaviors in complex environments.

Humans exhibit structured patterns of memory recall, including a tendency to recall more recent information and to recall events in the same order they were experienced. Classic computational models explain these patterns by positing that memories incorporate the ongoing ''temporal context'', formed by smoothly integrating the stimulus history. However, it is unclear whether a single mechanism can account for the full repertoire of human memory strategies, as the optimal approach may be task-dependent. For example, human memory experts widely apply the ''memory palace'' strategy, which is empirically better but not captured by temporal context models. Here we show that neural networks optimized for free recall develop diverse retrieval strategies, with only some of them resembling temporal context models.The best-performing models discovered a stimulus-invariant index code that emphasizes the studied position of each list item, instead of its temporal context. This creates a stable scaffold for forward recall akin to the memory palace technique. This index code was more likely to emerge when networks were i) encouraged to recall all studied items rather than prioritizing a few items, and ii) prevented from relying on recency, resonating with human data. Our findings demonstrate that human-like recall patterns can arise from multiple distinct computational mechanisms, and that sequential retrieval using item index is an optimal strategy that explains expert-level recall performance.

Building upon foundational psychological theories of event segmentation, this study addresses the limitation of overreliance on temporal boundaries as the primary segmentation criterion. Drawing on two experiments of direct and indirect causation in Mandarin Chinese, this study demonstrates how cognitive segmentation granularity and semantic integration jointly shape syntactic encoding. Results reveal distinct event encoding patterns for direct and indirect causation: coarse-grained segmentation leads to compact syntactic structures (e.g., verb-resultatives), while fine-grained segmentation yields varied multi-clausal expressions. Chinese speakers update event models via prediction errors of intentionality and protagonists, and tend to establish event boundaries at goal-relevant action endpoints when construing causal chains. These conceptual dimensions exert a modulating influence on both event segmentation and semantic integration. We propose a triad model integrating event segmentation, semantic integration, and linguistic specificity, providing a unified framework for elucidating the mind-language interface in conceptual construction and event coding of causation.

Classification of motor imagery (MI) tasks through EEG is valuable in brain-computer interfacing and rehabilitation engineering. EEG channels selection for MI task classification is well discussed problem and is challenging due to its combinatorial nature. Most of the existing methods are subject and task-dependent. This paper introduces a subject-independent EEG channel selection. The proposed approach consists of two stages. First, we rank channels based on their divergence from a reference channel Cz. We hypothesize that channels less divergent from Cz are more relevant for MI task classification. In the second stage, we employ a three-stage feature selection and classification model to evaluate the selected channels. It consists of a bandpass filter, followed by common spatial pattern (CSP) filter and three classifiers viz. SVM, 1-NN and 5-NN. Two publicly available datasets viz. PhysioNet and BCI Competition III IVa datasets have been used to assess the method. It performs 15.21\% more than 3Cs and just 2.91\% less than all-channels accuracy with as few as 20/118 channels on BCI Competition data and 19.64\% more than 3Cs on the PhysioNet dataset with 16/64 channels. Empirical comparison implies that the method performs better than classical models such as CSP Rank, fishers rank, and normalized mutual information, significantly. Results support that our hypothesis that divergence between channels and a reference channel Cz can be used as a ranking measure for channel selection.

The brain requires coordination among different regions to execute cognitive tasks, which may involve both positive- and negative-correlations. The topology of these correlations may indicate the mechanism underlying brain functioning in a given state. Here, we study changes in the functional connectomes (FCs) of both the positive and negative-correlations across various cognitive task states relative to the resting state, using publicly available electroencephalographic (EEG) data. Considering the EEG-specific topographical cortical regions as topographical modules (TMs), we find that the FC comprising positive correlations (G+) is modular. In contrast, networks of negative-correlations (G-) are anti-modular, with more connections between TMs than within them, and are associated with improved overall topological efficiency. These functional networks also show variability across frequency bands and brain states. In the low-frequency delta band, resting states exhibit higher modularity and anti-modularity than task states; in contrast, in the high-frequency Gamma band, modularity and anti-modularity are much higher during task states than in the resting state. The k-core analysis of all networks further reveals differences: G+ is more hierarchical and robust than G- across all states. Moreover, the task-state networks are always more hierarchical than the resting-state networks across all frequency bands. In the high-frequency gamma band, they are also significantly more robust than the resting-state networks. These networks also differ in the topology of their innermost core constituents: the innermost core regions of G+ are randomly connected and spatially localized, mostly in posterior brain regions across subjects, in the high-frequency gamma band. Whereas those in G- are spatially de-localized, cover the extreme anterior and extreme posterior brain regions, and remain anti-modular in all the frequency bands. Overall, our analysis reveals the presence of an anti-modular organization of functionally specialized TMs alongside their modular organization and points to task- and resting-state differences in their topologies.

Substantial progress has been made in recent years on understanding how the human brain parses and processes natural speech. Much of this progress has been based on modeling how brain activity relates to the different acoustic and linguistic features of speech. By fitting and testing models based on those features, one can test hypotheses about the kinds of computations and representations the brain uses to convert speech sounds into understanding. While much of this work has focused on modeling BOLD activity using functional neuroimaging or intracranially recorded electrophysiological signals, the approach has also proven useful with MEG and EEG. Indeed, noninvasive EEG has certain advantages for studying speech processing in terms of translational research and application. Research over the last decade or so has shown that EEG can be successfully modeled based on numerous acoustic, linguistic, and paralinguistic speech features. However, an important unanswered question hangs over all of this work: namely, what constitutes a good model of EEG responses to natural speech? Or, to put it another way, how much variance in EEG recorded during natural speech listening is explainable as having derived from that speech input? The present study aims to tackle this issue. We do so under the assumption that the best model for a person's EEG response to natural speech is a set of EEG responses from other people listening to the same speech. Using this assumption, we construct inter-subject models using EEG from 19 healthy adult native speakers of English who all listened to the same audiobook. The model for each subject involves predicting their EEG data using (dimensionality-reduced) EEG from different numbers of other subjects and then extrapolating to estimate the total explainable variance in the target individual's response to speech. Following this, we show that linear models (temporal response functions) based on several commonly used acoustic and linguistic speech features can predict most - but importantly not all - of the estimated total explainable variance in EEG responses across subjects.

Naturalistic cognition emerges from coordinated interactions among distributed brain systems operating across multiple representational scales. Characterizing this organization remains challenging because cognitively relevant information is embedded within high-dimensional neural activity. Here, we apply Multisubject Archetypal Analysis (MS-AA) [1] to naturalistic fMRI data collected during intact narrative listening, word-scrambled audio, and rest to investigate how condition-relevant information is distributed across archetypal representations. We examine both spatial and temporal formulations of MS-AA as complementary views of naturalistic brain activity. Across analyses, decoding performance consistently followed the hierarchy intact > word-scrambled > rest, indicating that archetypal representations preserve meaningful condition-related structure. Top-m decoding analyses further revealed that this information is highly compressible: relatively small subsets of archetypes frequently recovered substantial fractions of full-model decoding performance. Spatial AA exhibited a stable sparse-decoding regime that persisted across representational scales. Across a broad range of matched component ratios, approximately 5-15 archetypes consistently captured disproportionate amounts of condition-relevant information. These same sparse subsets also organized subjects into condition-aligned clusters more strongly than expected from random archetype subsets, with the strongest joint decoding-clustering effects occurring repeatedly within an intermediate representational regime (K {approx}50 - 88). Network over-representation analyses revealed that informative archetypes were not isolated canonical networks but distributed mixtures of interacting systems. Across the highest-performing decoding- clustering configurations, default mode and frontoparietal systems were consistently overrepresented relative to network size, whereas visual and limbic systems were underrepresented. Together, these findings suggest that the archetypal motifs most informative for distinguishing cognitive states are sparse, distributed subnetworks enriched for higher-order association systems. More broadly, the results demonstrate that MS-AA provides a useful framework for studying the compressibility, geometry, and multiscale organization of cognitive brain states.

Intense aesthetic experiences are among the most complex responses arising from the interaction of mind, brain, and context. Observations from fMRI suggest that when viewers feel highly moved by artworks, the underlying neural states differ from those accompanying less intense responses, particularly through recruitment of the DMN. Using electroencephalography and Bayesian category-specific cumulative link mixed models, we investigated whether such putative peak aesthetic responses exhibit threshold-specific neurodynamics rather than linear scaling with intensity. Twenty-two Chilean participants viewed 113 diverse local artworks whilst rating how moved they felt on a four-point scale. We analysed both canonical oscillatory power ({theta}-{gamma}) and aperiodic components (offset and exponent) during the contemplation window and the post-elicitor window. Threshold-specific effects were found: spectral features differentiated the highest rating category from moderate responses, rather than scaling uniformly across all intensity levels. During artwork visualisation, power in the {beta}1 and {beta}2 bands, as well as the interaction of {beta}1 with the 1/f exponent, predicted the transition to the most intense response; during the post-elicitor window, the aperiodic 1/f exponent predicted the transition from very low to higher-intensity responses. Modelling individual differences in spectral signatures (in the and {gamma} bands) credibly improved predictive performance (approximate leave-one-out cross-validation; elpd_loo), suggesting that neural variability reflects meaningful mechanistic heterogeneity in aesthetic processing rather than mere noise. These findings speak to a broader question, how the brain marks the intensity of conscious experience, and, more specifically, support the hypothesis that being intensely moved constitutes a qualitatively distinct neural state, characterised by specific configurations of oscillatory dynamics and cortical excitability that modulate the transition from low and moderate to peak engagement.

Intrinsic neural timescales represent the characteristic duration over which information is maintained in neuronal circuits. Evidence suggests that neural timescales vary systematically across the cortical hierarchy, with shorter timescales in primary sensory areas and longer timescales in higher-order association regions. In previous studies, hierarchy has been defined categorically, anatomically, or from the principal gradient of resting-state fMRI functional connectivity derived using diffusion map embedding (DME). Here, we assign hierarchical position to individual human intracranial electroencephalography (iEEG) recording sites by projecting their MNI coordinates onto this embedding space, derived from Human Connectome Project resting-state fMRI data. We estimated neural timescales from resting-state iEEG recordings in adult neurosurgical patients (n=46, 25 female) by extracting the aperiodic component of the local field potential power spectrum using spectral parameterization. Timescales increased monotonically with hierarchical position and associated with two region of interest (ROI)-level measures of network topology derived from DME of participants' iEEG functional connectivity: ROIs with stronger mean functional connectivity exhibited longer timescales, as did ROIs functioning as hubs, defined by proximity to the center of embedding space. Finally, timescales varied with sleep stage, with slowest values during NREM and fastest during wake and REM. The hierarchical gradient present during wake and N1 was no longer detected during REM, N2, and N3 sleep, driven by a selective increase in timescales at lower levels of the hierarchy. This work presents a novel metric of hierarchy that can be applied to iEEG data, establishes a direct link between neural timescales, cortical hierarchy, and network topology in human iEEG, and demonstrates that this hierarchical organization is dynamically modulated by brain state.